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We study the number of clusters in two-dimensional (2d) critical percolation, TVr, which intersect 
a given subset of bonds, F. In the simplest case, when T is a simple closed curve, Nr is related 
to the entanglement entropy of the critical diluted quantum Ising model, in which F represents the 
boundary between the subsystem and the environment. Due to corners in F there are universal 
logarithmic corrections to iVr, which are calculated in the continuum limit through conformal in- 
variance, making use of the Cardy-Peschel formula. The exact formulas are confirmed by large scale 
Monte Carlo simulations. These results are extended to anisotropic percolation where they confirm 
a result of discrete holomorphicity. 



I. INTRODUCTION 

In the percolation process sites or bonds of a regular 
lattice (or a graph) are independently open with a prob- 
ability, p and one is interested in the statistics of clusters 
of connected sites. In the thermodynamic limit for di- 
mensions d > 2 there is a critical point at p = p c , at 
which a second order phase transition takes place. Crit- 
ical percolation is scale invariant and its scaling limit is 
also believed to be conformally invariant. This has been 
rigorously proved for some lattices 1 -. Conformal invari- 
ance has been used to predict different relations about 
the order-parameter profiles, correlation functions, cross- 
ing probabilities, etc. of critical percolation. Some of 
the conformal results have been subsequently derived by 
rigorous mathematical methods, such as by Schramm- 
Loewner evolution^. 

In percolation the questions one usually asks con- 
cern the distribution of clusters, or the properties of the 
largest ones, but comparatively less attention is paid to 
the total number of clusters, which is expected to scale 
with the volume of the system under consideration. To 
the leading volume term there are corrections due to 
boundary effects, such as surfaces, edges and corners. 
The total number of clusters and its correction terms are 
generally non-universal, thus they depend on the type 
of percolation (bond or site) and the details of the lat- 
tice. The only exception could be the corner contribu- 
tion, which we are going to study in this paper. 

To be specific we consider bond percolation in a 2d 
square lattice and denote by T a subset of bonds. Ini- 
tially we take this to be a simple closed curve consisting 
of straight edges and corners and will generalize later. 
We are interested in the number of clusters, Ny, which 
intersect T, in the scaling limit when T is large but still 
much smaller than the total size of the system. This 
problem is closely related to the entanglement entropy of 
the bond-diluted quantum Ising model which is defined 
in the same lattice and T represents the interface, which 
separates a subsystem from the rest of the system. 



In the following we calculate Ny by different methods, 
both analytically and numerically. In the scaling limit 
we consider the random cluster representation of the Q- 
state Potts model and in the Q — > 1 limit we calculate 
N-p. In particular we obtain a logarithmically divergent 
contribution due to corners and its prefactor is calculated 
through conformal invariance. Using different geometries 
of T we have performed large scale Monte Carlo simula- 
tions, the results of which are confronted with the con- 
formal prediction. We have also considered anisotropic 
percolation, in which the probabilities (p x ,p y ) that hori- 
zontal and vertical bonds are open are different. 

II. POTTS MODEL REPRESENTATION 

Bond percolation is well known to be related to the 
Q —> 1 limit of the Q-state Potts model, defined on a 
lattice with sites i = 1,2, ... ,n and m nearest neighbor 
bonds. The partition function is given by 

Z(Q)=^I|exp(W s , fl ) , (1) 

s (ij) 

where the sum is over all spin configurations, s = 
(sx, S2j • • • > s n ), Si = 1,2, ... ,Q, the product runs over 
all nearest neighbor bonds, 5 SiySj is the Kronecker sym- 
bol and K is the reduced coupling, which is the ratio 
of the pair interaction and the temperature. Using the 
identity: exp (K5 Si , Sj ) = 1 + t^^,^ with p = 1 - e~ K , 
the sum of products in Z is written in terms of the so 
called Fortuin-Kasteleyn clusters^, denoted by F. In F 
the edge of the lattice i, j is occupied, if a factor j^6 SitSj 
is present and in any connected cluster the spins are in 
the same orientation. For a given element of F there 
are N to t(F) < n connected components and M(F) < m 
occupied bonds so that the partition function reads: 

z(Q) = J2 f Q Nm(f) (t^) 
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J2 F Q Nt ° t{F) P M{F) (l-py n ~ M(F) = (Q Ntot ) ■ (2) Usin^c=(6-/e)(3«-8)/2/sandv^=-2cas(47r/«) 



In this random cluster representation Q is a real param- 
eter and percolation is recovered in the Q — > 1 limit, 
when the critical point in the square lattice is given by 
p c = 1/2. The mean total number of clusters in percola- 
tion is: 



(N t< 



d]nZ(Q) 



dQ 



(3) 



Q=i 



Now impose a boundary condition that all the Potts spins 
are fixed (say, in the state 1) on T. The partition function 
is now: 



Zr(Q) ~ (Q 



Ntot-Nr 



>> 



(4) 



where iVr is the number of clusters which intersect T. 
Hence 



(N tot - N r ) 



ainZr(Q) 



dQ 



(5) 



3=i 



On the other hand at the critical point p = p c we can 
write: 



lnZ(Q) 
lnZr(Q) 



Af b (Q) 
Af b (Q) 



L T fs{Q) + Cv{Q)\nLT , (6) 



where A oc n is the total area, Lr is the length of T, 
and fb and f s are the bulk and surface free-energy densi- 
ties, respectively, which are non-universal. The last term 
represents the corner contribution. Hence we obtain: 



(N r ) = - f s (l)L r - C T (l)]nL r 



Cardy and Peschel^ considered these corner contributions 
to the free energy of a general conformally invariant sys- 
tem in domains with a boundary. Their results apply 
equally to an exterior boundary, with corners with inte- 
rior angle 7^, or an interior boundary, with 7^ replaced 
by 2tt — 7^. An important property of percolation is its 
locality, that is the partition function Zr(Q) is exactly 
the product of the partition function Z^(Q) for the in- 
terior of r, and Zp xt (Q) for the exterior. Thus we may 
apply the Cardy-Peschel result to each of these. Note 
that this would not be correct if we had given the clus- 
ters intersecting V a weight ^ 1. We therefore find from 
Ref. that the prefactor of the logarithm in Eq.© is 
given by: 



Cr(Q) 



24 ^ 

k 



27T - 7fc 



— ) _ (Ih 

Ik J ^ 7T 
27T - 7fc 
7T 



(8) 



where jk is the interior angle at each corner, and c(Q) 
is the central charge of the Q-state Potts model, and the 
two sets of terms come from the interior and exterior 
contribution^. 



we have: 



c'(l) 



5V3 

4-7T 



(9) 



The corner contribution to (iVp) which is derived here for 
bond percolation is expected to be universal, thus to be 
valid for site percolation, too. We are going to compare 
the conformal results in Eqs. (|7|8[ ) and ((9J with those of 
numerical calculations for different forms of T in SeclIVI 



III. ENTANGLEMENT ENTROPY OF THE 
DILUTED QUANTUM ISING MODEL 

The problem studied in Sec|H]for percolation is closely 
related to the entanglement entropy of a bond-diluted 
quantum Ising model and here we follow Refs. 0H . The 
model is defined by the Hamiltonian: 



(ij) ' 



(10) 



in terms of the a* ,z Pauli matrices at site i. The first sum 
in Eq.(|10p runs over nearest neighbors and the Jjj cou- 
pling equals J > with probability p and equals J = 
with probability 1 — p. At the percolation transition 
point, Pc, for small transverse field, h, there is a line 
of phase transition the critical properties of which are 
controlled by the percolation fixed points. The ground 
state of % is given by a set of ordered clusters, which 
are in the same form as for percolation. Now consider 
a subsystem, A the boundary of which is represented 
by r and calculate the entanglement entropy between 
the subsystem and the environment, which is given by 
St = —Tr (pa log 2 Pa) in terms of the reduced density 
matrix pa- Here we note, that in each ordered cluster 
the spins are in a maximally entangled GHZ state, thus 
for a given realization of disorder in the small h limit all 
those clusters give a unit (1) contribution to the entangle- 
ment entropy, which intersect T and contain also at least 
one point of the environment. We shall call them crossing 
clusters, thus 6>r is given by the number of crossing clus- 
ters and we are interested in their average value: (<Sr). 
Note that S r < N r , and thus (5 r ) < (N r ). The singu- 
lar corner contributions to (Nr) in Eq. ([7]), however, are 
due to large clusters, which are present typically among 
the crossing clusters, too. Consequently we expect the 
asymptotic form of the average entanglement entropy to 
be 



(Sr) = aLr + bhiLi 



(11) 



where a < —f' s (l) and b = — C T (1). The leading term 
in Eq.lfTTjl represents the so called area law, to which 
there is a logarithmic correction, which is expected to be 
universal and given in Eqs. (|7I8[) and ©■ 
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IV. NUMERICAL RESULTS 

We have performed large scale numerical calculations 
for site and bond percolations at the critical point on the 
square lattice 1 ^. The finite systems we have used have 
Lx L sites with L up to 8192 and with periodic boundary 
conditions. The different shapes of the subsystems used 
in the numerical calculations are illustrated in FigfT] In 
all cases the linear size of the subsystem is of O(L), so 
that Lr ~ L and in the following we shall use L instead 
of Lr- For a given shape of the subsystem and thus for 
its boundary T we have calculated the number of cross- 
ing clusters, Sr, which is then averaged: i) for a given 
percolation sample over the positions of T (typically 10 3 
positions), and ii) over different samples. Typically we 
have used 10 5 samples for each size, L, except the largest 
ones, where we had at least 10 4 samples. From the nu- 
merical data on (Sr) we have deduced estimates for the 
prefactor b by the so called difference approach. Here we 
use the relation: 

AS r (L) = 2 (Sr) (L) - (S r ) (2L) = blnL + const . (12) 

If r has a special shape, such as a square or a sheared 
square, (see the three shapes in the first row of Fig[lJ) 
we can calculate the corner contribution to (Sr) directly, 
by comparing results obtained at different geometries. 
This typ e of geometrical approach, which has been used in 
Ref. [tl|, will be explained in more detail in the following 
subsection. 




FIG. 1: (Color online) Shapes of the subsystems used in 
the numerical calculation: square, sheared square, anisotropic 
square; equilateral triangle, circle and section of a circle. In 
a lattice we consider the nearest neighbor bonds which are 
intersected by the contour of the subsystem and for each bond 
one of its sites belongs to F (depending on the local orientation 
of the surface) . 



A. Square subsystem 

The first geometry we consider for T is a square of lin- 
ear size I = L/2. The cluster structure of critical site 




FIG. 2: (Color online) Critical site percolation on a L = 128 
lattice. The subsystem is a square of size £ = L/2 = 64 and 
in this example iVp = 26 clusters intersect T = □ and the 
number of crossing clusters is Sq = 22. 

percolation is illustrated in Figj2] in which one can iden- 
tify the clusters which intersect T and also the crossing 
clusters. Here we use the geometrical approach, in which 
the complete square is divided either to four neighbor- 
ing squares with £ = L/2 or to orthogonal stripes of size 
L/2 x L: altogether four stripes due to the two different 
orientations. (This is illustrated in RefflTI in the right 
panel of Fig.l.) The two subsystems of different shapes 
have the same total boundary, however for stripes - due 
to periodic boundary conditions - there is no corner con- 
tribution. This is then obtained from the difference: 

Sg(L)=S n (L)-S str (L) . (13) 

Note that by this geometrical method the corner contri- 
bution is calculated for each sample, therefore the av- 
erage values have considerably less noise, than by the 
difference method. 

Making use of the fact that asymptotically (S^(L)) = 
61ni+const, we have calculated effective, size-dependent 
estimates for the prefactor by two-point fits, by com- 
paring the average corner contributions in finite systems 
of size L and 2L. The results are presented in FigJ3] 
both for site and bond percolations. With increasing 
L the effective prefactors approach a common, univer- 
sal limiting value of b = —0.077(1). This is to be com- 
pared with the conformal result in Eqs. (|7I8I) and © 
: -C"(l) = -573/(3671-) = -0.07657, thus the agree- 
ment is satisfactory. We note that a previous numeri- 
cal estimate in Ref. HI for smaller systems has obtained: 
b = —0.06(1). The corner contribution is related to the 
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FIG. 3: (Color online) Estimates for the prefactor of the 
corner contribution (Sfj) of the square subsystem obtained 
from two-point fits. The asymptotic values as L — > oo are the 
same for site and bond percolation and agree with the con- 
formal prediction. Inset: Large scale topology of the clusters. 
The square is divided into four quadrants and the different 
topologies of clusters in the quadrants are illustrated: empty 
box — > no site; filled box — > at least one occupied site. 



FIG. 4: (Color online) Ap = p — p c -dependence of the corner 
contribution to So for different sizes, L — 16, 32, . . . , 2048 
from below for bond percolation. Close to p c the extrapolated 
curve has a logarithmic singularity, see Eq.(fH]). In the inset 
the effective prefactors, b', are shown as a function of |Ap|. 



large scale topology of the clusters, as illustrated in the 
inset of Fig(3]and discussed in Sec lIVDI 

We have also studied the p-dependence of the corner 
contribution to («?□) outside the critical point by the ge- 
ometrical approach. As can be seen in Figl^S^" has a 
peak around p = p c and close to p c the extrapolated 
curve can be well described by the scaling result: 



Sq (p) ~ b' ln(p c — p) + const , 



(14) 



where b' = bis, with v = 4/3 being the correlation length 
critical exponent for percolation. Indeed, assuming the 
form in Eq. (fT4"|) we have calculated effective, p-dependent 
prefactors by two-point fits, which are shown in the inset 
of Figg] The extrapolated value for p -> p c (0.106(10)) 
is consistent with the scaling prediction (0.102). 



B. Sheared square subsystem 

The square subsystem used in the previous subsection 
is sheared now to a parallelogram, having the same sur- 
face: £ 2 = L 2 /4 and its (smaller) angle is 7. The geomet- 
rical approach to calculate the corner contribution to S^ r 
can be extended in this case for specific values of 7 given 
by the condition: tan(7) = 1/n, with n = 0, 1, . . . , 00 be- 
ing an integer. As shown in Figf5]with increasing L up to 
L = 4096 the numerical results approach the conformal 
prediction: 



C"(l) = - 



12 



I - 

,7 



1 



1 



1 



7T + 7 2TT 



7 
(15) 




FIG. 5: (Color online) Estimates for the prefactor —b for 
sheared squares as a function of 111(7). The numerical results 
are for L = 128, 256, . . . , 4096 from below, the full line 
represents the conformal result in Eq. (|15l) . In the inset the 
ratio of the extrapolated numerical data and the conformal 
result is given. 



Performing an extrapolation with an inL/L correction 
term we have obtained an excellent agreement, as shown 
in the inset of FigfSJ For other values of 7, which do not 
fit to the geometrical approach we have made calcula- 
tions by the difference method. Also in these cases the 
numerical results are found to agree with the conformal 
prediction. 
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C. Anisotropic percolation 



D. Corner probability 



Here we consider anisotropic bond percolation, in 
which the probabilities are p x (horizontally) and p y (ver- 
tically) and the critical point is given by the condition: 
Px + Py = 1- The system is a diagonally placed square 
with 2L 2 sites, and the subsystem is also a diagonally 
oriented square, having L 2 /2 sites and its boundary con- 
tains 4 x L/2 sites, see the third figure in the first row 
of FigJTJ This anisotropic system with symmetric angles, 
7r/2, is conjectured to be equivalent in the scaling limit to 
an isotropic system with asymmetric angles, 7 and 71 — 7, 
such that 



Py = sin((7r - 7)/3) 
p x sin(7/3) 



(16) 



This follows from the requirement that there exists a dis- 
cretely holomorphic observable^. More recently Grim- 
mett and Manolescui^ have proved that many properties 
of the scaling limit are the same if a more general inho- 
mogenously anisotropic lattice is embedded in the plane 
according to this prescription. 



For the (sheared) square subsystem we used the geo- 
metrical method, in which the L x L system is divided 
for four L/2 x L/2 squares and for four L/2 x L stripes 
and the difference in the number of crossing clusters in 
Eq. (|13[) is just the corner contribution. Having a general 
(connected) cluster it could occupy different parts of the 
four quadrants and its topology could be of four different 
types as illustrated in the inset of Fig(3J Among these 
the I, II and IV type of topology gives identical contribu- 
tion both for stripes and squares. Clusters with topology 
III, however, are crossing clusters for all the four possi- 
ble stripes, but these are crossing clusters only for three 
out of four squares. (In the left-bottom square there is 
no crossing.) Let us denote the occurrence probability 
of type III cluster as Ps(L), which is given by the ratio 
of such clusters in a L x L square, which have points in 
three quadrants, but have no point in the fourth one. As 
argued above Pz{L) is proportional to the corner contri- 
bution of crossing squares, more precisely: 



P 3 (L) = -5g"(L) = — -C(l)lnL + const. 



(17) 
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FIG. 6: (Color online) Estimates for the prefactor —b 
for anisotropic percolation for finite systems with L = 
32, 64, . . . , 1024 from below. The opening angle of the 
equivalent sheared square, 7 is given by Eq. (|16p and the full 
line represents the conformal result. In the inset the ratio of 
the finite-size results and the conformal results are given. 



This results is valid for sheared squares with angle 7, 
as well as for anisotropic percolation with a square sub- 
system. In these cases the appropriate results of S^'(L) 
have to be used. Interestingly this corner probability has 
a logarithmic L-depcndcnce and its prefactor is known 
exactly. 



E. Other subsystem geometries 

We have also studied subsystems with different geome- 
tries: equilateral triangle, circle and section of a circle, 
these arc illustrated in the second row of FigfT] In these 
cases the largest linear scale of the subsystem is fixed to 
L/2, while an angle was varied. In the calculations the 
difference approach in Eq. (|T2"j) has been used. 



1. Equilateral triangle 



In the anisotropic system we have chosen the prob- 
abilities in such a way, that the corresponding angle, 
7, has satisfied the condition: tan(7) = 1/n, with 
n = 0, 1, . . . , 00. In this way we could directly compare 
the results for anisotropic systems with those obtained for 
sheared squares in Sec. IIVBI The numerical results for 
the prefactor obtained on finite systems up to L = 1024 
are are shown in FigfB] For the three largest systems 
there is no systematic size-dependence of our data, the 
error is purely statistical and the numerical results agree 
well with the conformal prediction. 



For an equilateral triangle with a base angle 7 the ex- 
trapolated results for —b are shown in Fig[7] , which is 
compared with the conformal prediction: 



C(l) 



c'(l) 
24 



6 



1 



1 



7T — 27 



2 7 



2tt — 7 
(18) 

As seen in Fig. [7] there is a satisfactory agreement, al- 
though the statistical error of the numerical results is 
comparatively large, in particular for small and large an- 
gles. 
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FIG. 7: (Color online) Prefactor of the logarithm for an 
equilateral triangle as a function of the base angle 7 com- 
pared with the conformal result in Eq.lJTH]). The numerical 
results are extrapolated up to L = 4096 for site percolation. 
The statistical error is larger, than the difference between the 
conformal and the numerical results. 



2. Circle and section of a circle 




1 



7 < 7r there is a logarithmic correction, the prefactor of 
which is given by conformal invariance: 

12 [ \7 27I "-7/J 

The numerical results in Fig |5] are in agreement with the 
conformal prediction, although the statistical error of the 
numerical results is comparatively large. 

F. Line segment 

According to the derivation in Scc|Il] V docs not have 
to be a closed curve. For example it could be a straight 
line of length £. In that case we have only a corner con- 
tribution from two exterior angles, each 7 = 27r, so that 
Cr(Q) = -c(Q)/8. 

To study this problem we can use the geometrical ap- 
proach, when the line is oriented parallel with one of the 
axes of the square lattice. Then the corner contribution 
is related to the difference between the cluster numbers 
obtained for a periodic line of length L and that of two 
segments of lengths £ = L/2. In this case the corner con- 
tribution is simply half of the number of common clusters 
between the two line segments. 




4 5 6 7 8 9 10 
ln(L) 



FIG. 8: (Color online) Prefactor of the logarithm —b for a 
section of a circle as a function of the angle 7. The full line 
represents the conformal result in Eq. (|19[l , The differences 
between the conformal and numerical results are of the order 
of the statistical error. In the inset the corner contribution is 
shown for a circle: there is no systematic radius dependence. 
The numerical results are extrapolated up to L = 4096 (in 
the main figure) and up to L — 8192 (in the inset) for site 
percolation. 

According to the conformal results, there is no loga- 
rithmic correction to Sq for a circle shaped subsystem. 
Indeed the numerical results in the inset of FigJS] are in 
agreement with this statement: A5q approaches a finite 
limiting value of 0.65(5). 

In contrast, for a section of a circle with an angle < 



FIG. 9: (Color online) AS for the line segment obtained by 
the difference approach for different lengths, L. In the inset 
the p-dependence of S cr is shown close to the critical point 
for finite systems with L = 64, 128, 256 and 512 from below. 

In the numerical calculation we have used the differ- 
ence approach for a line segment of length I = L/2, which 
is placed in random positions and random orientations 
with respect to the axes of the periodic square lattice 
of size L x L. The average of AS\ (L) is shown in Fig. 
IH1 as a function of In L. The logarithmic dependence is 
clearly visible with a prefactor estimated as b = 0.086(1), 
which agrees fairly well with the conformal prediction: 

-Cri 1 ) = ~ 0.08615. We have also checked the 
p-dcpcndcncc of S? r (p), which is shown in the inset of 
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Fig. [9] A singularity is developed at the critical point as 
L — > oo. 

We note that if the line segment is put to the open 
boundary of the system, then the prefactor is differ- 
ent, being b = ps 0.1378 as obtained by conformal 
invariance&i^ and numerically b = 0.15(2)^. 

V. DISCUSSION 

In this paper we considered the logarithmic terms in 
the mean number Np of clusters in critical percolation 
which intersect a curve T in the cases when it has sharp 
corners or end points. We have shown that the Cardy- 
Peschel result^ can be simply applied to compute the uni- 
versal coefficients of these logarithmic terms, and that ac- 
curate numerical estimates agree very well with this, for 
a variety of shapes for T. We also considered anisotropic 
bond percolation on the square lattice and showed that 
if Eq. (TO)]) is used to deform the lattice then the predic- 
tions again agree with numerics if the correct effective 
corner angle is used. We have also pointed out a relation 
between the corner contribution to Nr and the statistics 
of cluster shapes. 

Our study is related to recent investigations of shape 
dependent terms of different thermodynamic quantities, 
mainly the free energy, of 2d critical systems. For crit- 
ical percolation, however, there is no logarithmic corner 
contribution to the free energy, since the central charge 



is c(Q = 1) = 0. In the latter problem the corner con- 
tribution to the cluster numbers and its higher moments 
are of interest, which are related to derivatives of c{Q) at 
Q = l. 

The excellent agreement of the numerical data with 
the theory serves as confirmation that critical percola- 
tion is indeed conformally invariant, of the Coulomb gas 
predictions for c'(Q = 1), and of the formula (fTo]) . Ex- 
tensions are obviously possible: for example by taking 
further derivatives with respect to Q we can find results 
for the corner contributions to the higher moments of 
the distribution of the cluster numbers Np. We note 
that recently^ an explicit example has been found of a 
correlation function in percolation which contains a mul- 
tiplicative logarithm. As with our result, this may be 
understood^ from the necessity of having to take a suit- 
able derivative with respect to Q at Q = 1 in the Potts 
model. 
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